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Abstract 

We compare the BFKL prediction for the associated production of forward jets at 
HERA with fixed-order matrix element calculations taking into account the kinematical 
cuts imposed by experimental conditions. Comparison with HI data of the 1993 run 
favours the BFKL prediction. As a further signal of BFKL dynamics, we propose to 
look for the azimuthal dependence of the forward jets. 



1. Introduction: The BFKL Pomeron has recently gained much interest, in particular in 
connection with the rise of F% at small Bjorken-x (xb) observed at HERA. Since it is generally 
believed that, at small x, F2 is driven by the gluon density in the proton, it is tempting to use the 
BFKL evolution equation for calculating the input distribution g(x, Qq). However, such quantita- 
tive BFKL-based predictions for xg(x, Q 2 ) are plagued by theoretical difficulties, and it is therefore 
important to look for other signals of the BFKL dynamics. A proposed clean "BFKL footprint" at 
HERA is the measurement of inclusive jets in deep inelastic scattering in the forward direction Q 
with longitudinal momentum close to the proton, i.e. the fractional momentum of the jet Xj e t 3> xb, 
and transverse momentum, kt, of the order of the virtual photon mass Q 2 . Analytic calculations 
[]|] and recent data from HI |3j show encouraging agreement, but the comparison still suffers from 
several shortcomings. First, the existing theoretical estimates of the cross section do not take into 
account the full experimental cuts, which are actually used in the data analysis. Secondly, in or- 
der to provide a solid "proof of existence" of the BFKL Pomeron, it is necessary to compare the 
data not only with the BFKL calculation but also with the fixed-order jet production cross section 
based upon QCD matrix elements. The BFKL prediction — which, as a result of gluon production 
between the forward jet and the current jet, rises as a function of 1/xb — is expected to lie above 
the matrix element calculation, and the data will discriminate between the two predictions. A third 
obstacle in verifying the BFKL prediction is the lack of hadronization in the theoretical analysis: 
Monte Carlo studies indicate that, for this type of measurements, hadronization may not be too 
strong an effect, but there does remain a serious uncertainty. An improvement in this direction 
might be achieved by the development of a BFKL-based Monte Carlo program. Further evidence 
for the BFKL Pomeron can be obtained by looking into the azimuthal dependence of the forward 
jet cross section: in the HERA reference frame BFKL dynamics predicts Q that, with increasing 
rapidity difference between the forward jet and the current jet, the forward jet "forgets" about the 
azimuthal direction defined by the outgoing electron, and the cross section becomes ^-independent. 
This decorrelation is, again, a result of the radiation of gluons between the forward jet and the 
current jet. Correspondingly, in the fixed-order matrix element calculation one expects the <3?- 
dependence to persist, even for large rapidity intervals. For the Mueller-Navelet Q jet production 
at hadron colliders this decorrelation effect has been investigated JjJ , and it has been demonstrated 
that such a measurement might give a clear signal for the BFKL Pomeron between the two jets. 
For the analysis of experimental data from HERA an analytic calculation was performed recently 
[]|]. A numerical analysis, however, has not been done yet. 

In the present paper we will first improve on the comparison between the cross section formula 
and the data. Based upon the analytic expressions given in Ref. H we compute the cross section, 
taking into account the full experimental cuts used in the recent HI analysis Q. We then make 
use of the fixed-order matrix element calculations for the two- and three-jet cross section || and 
compute, with the same phase space cuts, the numerical predictions. Both calculations are then 
compared with the data. In the final part of this paper we suggest a new measurement of the 
azimuthal decorrelation to reveal the BFKL Pomeron. We discuss and present a numerical anal- 
ysis of the decorrelation in the azimuthal angle <£, both in the BFKL framework and within the 
fixed-order jet cross section calculation. 

2. Experimental Considerations: The study of associated jet production is an experimen- 
tal challenge. It was shown in Ref. || that the requirement of Xj e t/xs to be large results in typical 
jet angles of a few degrees with respect to the proton direction, also termed forward direction. Due 
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to the unavoidable beam-pipe hole for the detectors at the ep collider HERA, the acceptance is 
limited to jets with an angle larger than, for example, 6° in the present HI detector. At smaller 
angles the jets are insufficiently contained in the detector and the experimental separation from 
proton remnant fragments is difficult. The angular requirement leads to maximally reachable x jet- 
values of about 0.1. Furthermore, due to the limited amount of data, a compromise has to be found 
for the condition k 2 /Q 2 ~ 1. This condition is essential since it defines a kinematical environment 
where no ambiguity between the DGLAP evolution and the BFKL Pomeron exists. In practice an 
interval is defined around this central value of k 2 /Q 2 . The HI collaboration has made a first mea- 
surement of the associated or forward-jet production in deep inelastic scattering Ml. The selections 
are guided by the limited statistics (300 nb _1 ) and can be refined for future large-statistics samples. 
In this paper we adopted similar selection criteria. Deep inelastic events are selected by requiring 
the energy and angle of the scattered electron to be larger than 12 GeV and less than 173°, respec- 
tively. The first requirement avoids regions with large contributions of radiative DIS events and 
experimental background due to misidentified photoproduction events, while the second matches 
the acceptance of the calorimeter that detects the scattered electron. In order not to confuse the 
forward jet with the one at the top of the ladder, the requirement y > 0.1 was imposed to ensure 
that the jet of the struck quark is well within the central region of the detector and is expected to 
have a jet angle larger than 60°. Experimentally a cone algorithm is used to find jets, requiring an 
E t larger than 5 GeV in a cone of radius R = y/ Arj 2 + Acft 2 = 1.0 in the space of pseudo-rapidity r\ 
and azimuthal angle in the HERA frame of reference. For the jet selection we consider two sets of 
cuts, the first set (selection A) matches the published data, while the second set (selection B) could 
be used when more statistics becomes available and the systematics gets under better control. For 
selection A jets are accepted as forward jets if Xj e t > 0.025, 0.5 < p 2 et /Q 2 < 4, 6° < 6j e t < 20° 
and pj e t > 5 GeV, where 9j e t is the forward jet angle, pj e t is the transverse momentum of the 
jet, in approximation of k t . Selection B accepts forward jets if Xj e t > 0.035, 0.5 < p 2 et /Q 2 < 2, 
7° < 6 j e t < 20° and pj e t > 3.5 GeV. These selection criteria allow a study of the cross section of 
forward jet production in the region Q 2 « 20 GeV 2 and 1 x 10~ 4 <i B <2x 10~ 3 . Hence the ratio 
Xj e t/xB is always larger than 10. 

3. The BFKL Approach: The BFKL approach makes a prediction for the high-energy (xb <C 
Xj e t) behaviour of the cross section for the process g — > g + (ng) + qq or q — > q + (ng) + qq. 
Consequently, the first nonvanishing contributions appear in the (3+l)-jet matrix elements^] aris- 
ing from graphs with one gluon in the i-channel (for the (2+l)-jet processes q — > qg and g — > qq, 
there are no contributions with gluons in the i-channel, and therefore their dependence on Xj e t/xB 
is power-suppressed). A generic diagram is shown in Fig. la. 

Using the gauge condition (q + xp) ■ A = 0, where q and p are the four-momenta of the virtual 
photon and proton, respectively, the diagrams corresponding to the two processes g — > gqq and 
q — ► qqq have to be calculated. The incoming parton is always scattered into the forward region, 
well separated in rapidity from the photon region. Moreover, the two diagrams differ only in their 
colour content, since the coupling of the i-channel gluon to the incoming parton (gluon or quark) 
is of eikonal type. Due to the large rapidity gap given by ln(x jet/xs) the probability to emit more 
gluons is increased, for the smallness of the strong coupling constant a s is compensated by the 
large logarithm \n.(xj e t/ xb)- This means that higher-order corrections become relevant, and in the 
end all diagrams with an arbitrary number of s-channel gluons have to be considered. A generic 

1 The notation "(n+1)" stands for n jets in the current fragmentation region and one jet in the target fragmentation 
region, the latter consisting of the proton remnants. 
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Figure 1: Generic diagrams for the process (a) g — > g + qq and (b) g — > g + (ng) + qq. Diagrams 
with an incident quark (not shown) have been considered as well. 



diagram of this process is shown in Fig. lb. Their resummation is performed in terms of the BFKL 
equation [[jj and leads to the following result for the one-parton inclusive cross section |J: 



with 



\g(x jet ,k 2 ) + ^ ^q f {x jeU k 2 t )) (1) 



w/(») \^ n2 njFa 11 f°° a I + ^ sinh(//7r) 



co T ln( g 



co S (mMt4) » ( 2 ) 



XB wl n) = Y.Q}*s s[m> \ /f d/M \ 



00 j + fJ 2 sinh(//7r) 



1 + /x 2 /i cosh (/j,7t) 

. (^)" M cos (, M |)) (3 ) 

where 

W (n )M ) = ^{^(1) - ^(^ + ^) " ^( ! 4 Ll -^)} • W 

Here n > is an integer, and Wt and Wi refer to the transverse and the longitudinal parts of the 
hadronic tensor W pa . The transverse momentum of the jet and its longitudinal energy fraction are 
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given by kt and Xj e t- The first-order running coupling constant was used with kt as scale. The main 
features of the BFKL prediction, in particular the power-like increase at small xb, can immediately 
be read off from Eqs. (Q) and (|3|), making use of the fact that the integrand dominates at /x = 0. For 
/i = and n = one finds 9/2 for the ratio of the transverse and the longitudinal part of the cross 
section. For large \xi{xj e t/xB) the BFKL power is obtained from u(0,0). Before we discuss our 
numerical results based upon this BFKL formula, we briefly turn to the fixed-order matrix elements. 

4. Fixed-Order Matrix Elements Calculations: In order to compare the BFKL results with 
a fixed-order calculation based on exact matrix elements, we study the contributions for the pro- 
duction of two and three partons on the Born level, giving rise to (2+1)- and (3+l)-jet events. We 
use the PRO JET Monte Carlo program || , where the cross sections from Ref. |8| are implemented 
in the modified JADE jet-definition scheme^]. Two partons can be produced in the parton-level 
processes q — > qg and g — > qq; for three partons there are the processes q — > qgg, g — > qqg and 
q — > qqq. These matrix elements suffer from infrared and collinear singularities, which are excluded 
by means of suitable cuts on the parton momenta. This is done by requiring that the invariant mass 
squared s a b = ^PaPb be larger than some mass scale M 2 ut = cW 2 , c being the jet cut parameter and 
W 2 the squared invariant mass of the hadronic final state; a and b are set to q, q, g and r, where 
p r is the momentum of the proton remnant jet. All matrix elements except those from the class 
g — > qqg are singular if any of the invariants goes to zero. The matrix element g — > qqg is not sin- 
gular for s q q = 0, because of the additional propagator of the quark line. In the following numerical 
study, we have used the same cuts for the forward jet as those that are applied by the HI collab- 
oration in their BFKL study. The parton density used is the CTEQ3 leading-order parametrization. 

5. Numerical Results: Numerical results for both approaches are shown in Fig. 2. Let us 
begin with the BFKL approach. The plots depicted in Fig. 2a were evaluated by integrating the 
differential cross section within the given cuts (see above). The forward jet selection A is used. 
Two cases are considered: 



(i) The fixed-order but high energy asymptotic expression for the two processes q — > q + qq and 
g — > g + qq (lower two curves in Fig. 2a). These cross sections are evaluated in the BFKL-type 
high energy approximation, and from Eqs. (]2|) and (|3|) they are obtained by simply setting 
uj(n, /i) = 0. 

(ii) The full BFKL result according to Eqs. (g) and (||) (the upper two curves in Fig. 2a). Only 
light quarks (u, d, and s) are included. The curves were scaled by 20% to approximately 
include the effect of charm production, in accord with recent data |0|. 

Case (i) shows good agreement with the results of the second approach, the full matrix element 
calculations (Fig. 2b, see below). The full BFKL approach, however, scales up the normalization 
roughly by a factor (xjet/xs) ' 5 and reveals a much steeper shape in the x^-distribution. In 
comparison with data, it is in particular the shape of the cross section which is an important test 
for the presence of BFKL. The absolute normalization suffers from some uncertainty as discussed 
below. Figure 2c shows the same calculation for the second set of kinematical cuts. 

2 Another calculation of (3+l)-jet cross sections can be found in Ref. |E3]. It should be noted that the JADE 
jet-definition scheme is different from the cone scheme employed in the experimental analysis ym 
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Figure 2: (a) xb dependence of the cross sections from the BFKL calculation (incoming quark and 

gluon ( ) and incoming gluon only ( )) and from an approximate analytical calculation 

of the three-parton matrix elements (incoming quark and gluon ( ) and incoming gluon only 

(■■■))■ (b) xb dependence of the cross sections from the exact fixed-order matrix elements, with 
any of the partons as the forward jet. g — > qqg integrated over the full range of the invariant 

s q q ( ); the same, but with the additional invariant mass cut s q q > M^ ut ( ); q — > qgg 

(— ■ — ) and q — > qqq (— ■■ ), with the cut sab > ^cut on a ^ invariants; sum of q — > qg and g — > qq 

( ■■), with the cut sab > M^ ut on all invariants. The jet cut parameter is c = 0.005. The 

distributions (c) and (d) are the same as (a) and (b), but for the jet cuts set B and calculated for 
10~ 4 < x < 4 x 10- 3 . 
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Figure 3: Comparison of the BFKL calculation and the approximate analytical calculation of the 
three-parton matrix elements, shown in Fig. 2a, compared with the data of the HI experiment ^jj. 

Next we turn to the second approach, the fixed-order matrix element calculation. The results 
for the x^-dependence are shown in Figs. 2b and 2d for the jet selections A and B, respectively. As 
a result of the applied cut, the contribution from the (2+l)-jet matrix element is small. The (3+1)- 
jet contribution q — > qgg is small as well; the quark densities are small compared with the gluon 
densities, and the process is not a £-channel process. The q — > qqq contribution is larger, because 
in this class of diagrams there is a t-channel gluon which enhances the cross section for forward-jet 
production. The (3+l)-jet cross section is dominated by the process g — > gqq. The gluon density 
is large, and there are large contributions from the diagrams with a t-channel gluon. In order to 
compare the full fixed-order matrix element to the BFKL prediction, it is, in principle, necessary 
to calculate a cross section that is inclusive in all particle momenta except for the momentum of 
the forward jet. In the present version of PROJET, such a study is not feasible. However, it is 
possible to remove the cut s q q > M^ ut for the process g — ► gqq (see the above discussion). Relaxing 
this cut increases the cross section considerably. The shape of the xb distribution compares well 
with that of the approximate analytic calculation for the process g — > gqq (see above), whereas for 
the full BFKL result (Fig. 2a) the slope in is somewhat steeper. It should be noted that the 
cross section from the exact matrix elements are calculated from a modified JADE scheme and are 
less "inclusive" than those from the approximate analytic calculation; a direct comparison of the 
magnitude is therefore not possible. 

Finally, we compare the analytic predictions with the HI data. Because of the JADE jet definition 
scheme, the fixed-order matrix elements are not quantitatively comparable with the data based on 
the cone scheme. Moreover PROJET treats the production of charm quarks in the massless approx- 
imation. The analysis presented in Ref. § has used only two a^-bins, and the two values are shown 
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in Fig. 3. Obviously, the data are higher than the analytic fixed-order results, indicating that gluon 
production, which is responsible for the rise of the BFKL Pomeron, is present. On the other hand, 
the data points lie very close to the BFKL prediction. This looks very encouraging and could be 
interpreted as the first sign of the presence of the BFKL Pomeron at HERA. However, we feel that 
several remarks are in order. First, for the analogous process at hadron colliders (Mueller-Navelet 
jets) a comparison of a next-to-leading order QCD matrixelement calculation with the BFKL ap- 
proximation has shown [12| that the latter one overestimates the available phase space. A more 
accurate treatment of the forward jet production, therefore, would most likely tend to lower the 
BFKL prediction. A similar effect has to be expected also in the forward-jet production at HERA. 
Secondly, our treatment of the charm contribution is clearly unsatisfactory, and a more accurate 
treatment of the charm mass has to be included in our calculation. Finally, as we have mentioned 
in the beginning, also hadronization which is not included in any of our calculations may have some 
effect. We therefore feel that we have to be cautious in drawing too strong conclusions. It seems 
safe to say that - within the accuracy to which the data can be compared with our calculations - 
the measurement of forward jets does not agree with the approximate fixed-order matrix element 
calculation, and the comparison with the BFKL prediction is very encouraging. However, there 
may be some effects, which can distort this good agreement. 



6. Angular decorrelation: Recently || the azimuthal dependence of the jet production cross 
section has been calculated^]. The angle $ is defined in the HERA frame as the angle between 
the transverse momenta of the outgoing lepton and the forward jet, and the formula for the cross 
section reads: 



- 2(1 - y) cos(2$) x B wj 2) } (5) 



g(x jet ,kf) + - J2<lf( x jet,kt 



2^kf 1 9 



Wj is given by Eq. (|3|) with n = 2. For small 1/xb, Eq. dBI) predicts a maximum of the cross 
section at $ = tt/2. When 1/xb gets large, we expect in the BFKL prediction the ^-dependence 
to disappear. In Eq. (||) this effect follows from the difference of the powers u)(Q, 0) ~ 0.50 and 
cj(2,0) ~ —0.22. With increasing Xj e t/xB, the angular dependent part is suppressed whereas all 
other parts increase. 



Numerical results are shown in Fig. 4. In the BFKL case (Fig. 4a) the ^-dependence is completely 
washed out when x b becomes small (horizontal line in Fig. 4a) . This is in complete agreement with 
the theoretical expectation. For the approximate fixed-order process g — > gqq, on the other hand, 
the cross section is not flat. It has a clear maximum at <3? = tt/2 and is found to be rather insensitive 
to the choice of the xb bin (not shown). Figure 4c shows the same calculations as Fig. 4a, for the 
jet selection set B. 

For comparison, we also show the results of the fixed-order matrix element calculation (Fig. 4b). 
3 A somewhat different type of azimuthal dependence has been studied in Ref. [H. 



7 



1 .50 



.25 



1 .00 



0.75 



0.50 





(a) 




I 



0.00 3.14 6.28 



1 .50 



1 .25 



1 .00 



0.75 



0.50 




0.00 3.14 6.2^ 



1 .50 



1 .25 



.00 



0.75 



0.50 




0.00 3.14 6.21 



1 .50 



1 .25 



.00 



0.75 



0.50 





(d) . 




'14 



0.00 3.14 6.2^ 



Figure 4: Dependence on the difference <3? of the azimuthal angles of the outgoing lepton and the 
forward jet, normalized to a common average, (a) BFKL result, the bin in xb is [0.2 x 10~ 3 ,0.4 x 

1CT 3 ] ( ). Approximative analytical Born result, the bin in xb is [0.2 x 10~ 3 ,0.4 x 10~ 3 ] 

( ). (b) Fixed-order result, the bins in xb are given by [0.2 x 10 _3 ,0.4 x 10 -3 ] ( ), 

[1.8 x 10~ 3 ,2.0 x 10~ 3 ] (■■■jand[3.4xl0- 3 ,3.6xl0~ 3 ] ( ). The jet cut parameter c = 0.005. 

The matrix element is g — > qqg, integrated over the full range of the invariant s q q, where any of the 
partons may be the forward jet. The distributions (c) and (d) are the same as (a) and (b), but for 
the jet cuts set B. 
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We have used the jet cut M 2 ut = cW 2 with c = 0.005. The curve shown in Fig. 4b belongs to 
the process g — ► gqq, integrated over the full range of the invariant s q q. The forward jet may 
belong to any of the partons. At smaller xb, there is a tendency to enhance the maximum of the 
distribution at <E> = n/2 and $ = 3tt/2, similar to the BFKL Born calculation in Fig. 4a. There is 
a clear difference with the full BFKL prediction in Fig. 4a. At larger xb, on the other hand, the 
full matrix element calculation does not agree very well with the BFKL Born calculation. Figures 
4c and 4d show the results for selection B of the jet cuts. 

The present detectors at HERA allow for jet angular measurements with a resolution of about 
5°-10° for jets with an Ef larger than 5 GeV. This opens a window for measuring angular cor- 
relations between struck quark and forward jet, or between scattered electron and forward jets. 
The major background for this measurement is made up of radiative events; due to the changed 
kinematics at the hadron vertex, these can eject a jet in the forward direction. This background is 
however expected to be controlled experimentally to a sufficient level to allow for a measurement of 
better than a 10% significance level in large data samples. We therefore expect that a measurement 
of the angular decorrelation and the comparison with theoretical expectations provides a further 
test of BFKL dynamics. 

7. Conclusions: In this paper we have made a comparison for the production cross section of 
associated jets in the forward direction of the BFKL prediction and of the fixed-order matrix ele- 
ment calculation. Both curves show a clear difference at small xb, thus confirming the expectation 
that this process provides a clean signal of the BFKL Pomeron. Confrontation with the data is 
encouraging: the data clearly lie above the fixed-order calculation and prove that the radiation of 
extra gluons is an important effect. The BFKL calculation is very close to the data, which could be 
the first sign of the presence of the BFKL Pomeron in the HERA data. However, there are several 
aspects that need further study, in particular the contribution of charm and hadronization effects in 
the final state. From the theoretical side, there is also the question of the kinematic region in which 
the BFKL approximation is applicable. We expect that the strongest correction will come from 
unitarization effects, and we know that first unitarity corrections will lower the BFKL prediction. 
Hence it may be premature to interpret the data as a proof of the presence of the BFKL Pomeron. 

For future measurements we have suggested to look also into the azimuthal dependence of the 
jet cross section. Comparison with fixed-order calculations show that this provides a new signal for 
BFKL dynamics, which may help to find further footprints of the BFKL Pomeron at HERA. 
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